Sublattice model of atomic scale pairing inhomogeneity in a superconductor 
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We study a toy model for a superconductor on a bipartite lattice, where intrinsic pairing inho- 
mogeneity is produced by two different coupling constants on the sublattices. The simplicity of 
the model allows for analytic solutions and tests of the consequences of atomic-scale variations in 
pairing interactions which have been considered recently in the cuprates. We present results for the 
transition temperature, density of states, and thermodynamics of the system over a phase diagram 
in the plane of two pairing coupling constants. For coupling constants of alternating sign, a gapless 
superconducting state is stable. Inhomogeneity is generally found to enhance the critical tempera- 
ture, and at the same time the superfluid density is remarkably robust: at T = 0, it is suppressed 
only in the gapless phase. 

PACS numbers: 74.81.-g,74.62.-c,74.78.Na,74.25.Dw 



I. INTRODUCTION 



In recent years, the effect of underlying inhomo- 
geneities in superconductors has attracted many re- 
searchers. In the context of high temperature super- 
conductors (HTSC), checkerboard local density of states 
(LDOS) oscillations and strong nanoscale gap inhomo- 
geneity have been observed in scanning tunnelling spec- 
troscopy (STS) experiments 1 , and signals in dynami- 
cal susceptibility measured by neutron scattering have 
been interpreted as stripelike nanoscale modulations of 
charge and spin degrees of freedom 2 . While it is still not 
clear whether these modulations are intrinsic in nature or 
driven entirely by disorder, others have pointed out that 
the very existence of inhomogeneity may enhance T c , and 
implicitly suggested that the phenomenon of high criti- 
cal temperatures may rely on it 3,4 . In particular, Martin 
et al. studied a simple model, in which they included a 
ID inhomogeneous pair potential varying on an arbitrary 
length scale £ > a, where a is the lattice spacing, and 
argued for an increase of the critical temperature with 
inhomogeneity, with a maximum enhancement obtained 
for I ~ £oj where £o is the coherence length of the analog 
homogeneous system. Similar results were obtained by 
Loh and Carlson, where they found however that an in- 
crease in the transition temperature occurred at the cost 
of reduced superfluid density 5 . Aryanpour et al. 6 per- 
formed a systematic numerical study of s-wave supercon- 
ductivity on checkerboard, stripe and random patterns of 
pairing potentials in an attractive Hubbard model; they 
also found a positive correlation between critical tem- 
perature and inhomogeneity for many situations. It is 
also known that models with enhanced T c and signifi- 
cant inhomogeneity are associated with unusual features 
in electronic structure. For example, recent calculations 
of high-quality ultrathin s-wave superconducting films 
have shown that the surface-driven inhomogeneity can 
enhance T c due to weakly dispersive low-energy subgap 



states in the system 7 . The above studies were primarily 
considerations of BCS type models. Recently, however, 
Tsai et al. 8 have studied small size repulsive Hubbard 
models with a "checkerboard" -like modulation of hop- 
ping matrix elements using exact diagonalization, and 
also concluded that the inhomogeneous models could lead 
to pairing at very high binding energy scales. 

Nunner et al. 9 underlined the importance of this phe- 
nomenon for cuprates by explaining many of the correla- 
tions between the measured local gap and other obscrv- 
ables in Bi-2212 by assuming a <i-wave pairing interaction 
g(r) varying on a unit cell size, and driven by dopant 
disorder. If the effective pair interaction can vary on 
an atomic scale, and the BCS coherence length is suf- 
ficiently short, one can imagine situations where in a 
complex periodic crystal with large unit cell the pairing 
interaction varies over an atomic length scale in the ab- 
sence of disorder, i.e. a periodic, modulated g(r), which 
will give rise to a periodic, modulated order parameter 
A(r). The properties of such a system are important 
to understand, because it represents the simplest type of 
nanoscale pairing modulation, and one may study the re- 
lationship between superconducting properties and inho- 
mogeneity easily. Aryanpour et al. 6 considered systems 
of this type within a purely numerical approach, and ob- 
tained a variety of interesting results. Nevertheless, in 
such a treatment it is sometimes difficult to extract sim- 
ple conclusions and understand exactly which aspects of 
superconductivity are really enhanced or suppressed by 
inhomogeneity, and why. 

In this paper, we consider the simplest such model, 
a bipartite lattice in two dimensions with two different 
values of g on two interpenetrating sublattices. The ho- 
mogeneous version of the model has a superconducting 
ground state on a bipartite lattice in arbitrary dimen- 
sion and filling 10 . Here we calculate the ground state and 
thermodynamics of this system for any values (gA,gB), 
where A and B represent the two sublattices. We distin- 
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FIG. 1: T = mean field phase diagram in the space of two 
coupling constants §a and found in this paper. The solid 
lines represent transitions between normal metal (defined by 
= A_b = 0), gapped and gapless superconducting phases. 
The dashed-dotted line represents the line of constant average 
pairing g = t studied in the text. The dashed line represents 
the homogeneous BCS case with g A = g B . Finally, the dot- 
ted line represents the normal-gapless superconducting transi- 
tion line as calculated analytically within the approximation 
described in the text using the "window" density of states 
(wDOS). 



guish however two cases: one, where both g's have the 
same sign; and a second, where they differ. Situations 
like the latter case where the pairing potential changes 
sign abruptly over a unit cell distance have been con- 
sidered in the impurity context 11,12 , but not in periodic 
cases, and it is not clear whether true superconducting 
solutions can exist in such a situation. 

In section II, we introduce the sublattice pairing 
model and give exact expressions for the order param- 
eter and critical temperature T c within the generalized 
BCS scheme. We show that T c can indeed be enhanced 
by detuning the coupling constants from the BCS point 
9 A = 9b- We then derive in section III the quasiparticle 
excitations and thermodynamic properties of the model. 
When the average pairing strength is held constant but 
inhomogeneity is increased, the model displays a shrink- 
ing spectral gap in the total DOS. In the case of opposite 
sign coupling constants, the excitation spectrum becomes 
gapless, and the mean field solution corresponds to an 
order parameter which oscillates in sign from point to 
point. In section IV, we calculate the supcrfluid density 
n s for the model, and show that for the case of same-sign 
coupling constants inhomogeneity preserves a robust n s 
at T = 0. In the case of opposite signs, however, n s is 
suppressed by the onset of a residual normal fluid. The 
variations of superconducting phenomenology found in 
our solution to the model is summarized in the phase 
diagram shown in Fig. 1. 




FIG. 2: Unit cell of square lattice, where squares and circles 
denote sites of type "y4" and "B", and the lattice spacing is 
denoted by a, which we have set to 1 throughout our discus- 
sion. The coordinate system used in Fourier transformation 
is shown with two dashed lines. 



II. MODEL 

The Hamiltonian is given by Eq. (1), where Cj, c\ are 
the electron annihilation and creation operators, t is the 
strength of hopping energy between the nearest neighbors 
and A, is the mean field gap at site i. The sum in the 
Hamiltonian is over all sites and spins denoted by a. The 
unit cell of the bipartite lattice is shown in Fig. 2, and the 
Fourier transformation in Eq. (2) is defined with respect 
to this cell. We first consider the case of half-filling, for 
which the chemical potential /j, is set to zero. 

In anticipation of the doubling of the unit cell due to 
an interaction specified on each sublattice, we Fourier 
transform fcrmionic operators as 

c%,a = 51 c ki<T e lk ' r \ (2) 

k 

where the sum now runs over the arising for staggered 
ordering ( "antiferromagnetic" ) half-Brillouin zone speci- 
fied by < \ki\ < Tr/(y/2a) for i = x,y, with directions 
specified in Fig. 2. The Hamiltonian now becomes 

h = E^+A A 4M A , + A B ^4 B ff (3) 

k,<T 

= 5>t-A/-c k , (4) 

k 

where we have assumed s-wave pairing and the dispersion 
relation is given by 

£ fc = -4tcos(^)cos(% (5) 
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The matrix M may be expressed as 



M 



£fc A A 

Cfe A B 

A A 

A s 



(6) 



in the staggered particle-hole space spanned by Ck = 
(c^ k(T , c^ kcr , c^l^, c^l^). After performing the canonical 
transformation 



4:=< Q (khL+<<(k)7k,-„ 



(7) 



with a — A, B, we find the quasiparticle energies ±Ei^, 
with 



E 



TA A ± A 1 



{A A + A B ) A +^l 



1,2 



(8) 



which reduce to the usual E^ = \J £^ + A 2 for A" — > A. 

The eigenstates of the problem, of the form 
\u A , u B , v A , v B ), are given by 



\\x,~Ei >= [-xi,-x 2 ,-x 1 ,x 2 \ 

\\ 2 ,-E 2 >= [-X2,-X!,X2,-Xi] 
\h,E 2 >= [-X 2 ,Xi,-X2,-Xi] 

|A 4 , Ei >= [-xi,X2,x 1 ,x 2 ] 
where 



v 72 (*? + *1) 
i 

^2 (xl + 4Y 
i 

yj2(x( + xl)' 
1 



si,a = y y (A A + A B ) Z + 4£ 2 t (A a + A- 8 ). 
Using the definition of the gap functions, 



Af = g a = E^«<« th 



/ /3-En 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 



we can write down the gap equation for the A Q 



a a = j: 

k 



g- 



-x\ th 



13E X 



xl th 



x 2 th 



V 2 



xf th 



2 



(15) 



Note that within our convention a positive g corre- 
sponds to an attractive on-site interaction. We now cal- 
culate T c by linearizing the gap equations. To obtain 
analytical results, we estimate the integrals involved us- 
ing the approximate window density of states (wDOS) 




-At<uj < At 
elsewhere , 



(16) 




FIG. 3: Critical temperature T c /t plotted vs. difference of 



sublattice coupling constants S 



g B )/t for fixed g = t. 



Squares: window density of states. Circles: exact. Insert: 
T c /(4t) vs. g B /t for fixed g A = t. 



which is a good approximation for the tight binding 
model for qualitative purposes. Then T c takes the form 



k B T c 

9eff 



exp 



1 



9eff 



Ag A g B - 8i (g< 



8t(g A + g B )~(8ty 



(17) 



The dimcnsionless effective interaction <? e // plays the 
same role as the single attractive coupling constant in 
BCS weak-coupling theory. The phase transition to the 
normal state takes place as g e ff — * + , determining the 
dotted line in Fig. 1. Note that the mathematical ex- 
pression for T c (17) actually diverges as g e ff ~ * _ , i.e. 
as one approaches the phase transition from the normal 
state side, again in perfect analogy with BCS. The appar- 
ent nonzero value for T c on this side is of course spurious, 
as there is no ordering below this temperature, as will be 
shown below. 

This estimate for T c is compared in Figs. 1 and 3 with 
the exact numerical solution for the normal state insta- 
bility temperature from Eqs. 15. In order to judge the 
extent to which inhomogeneity enhances superconductiv- 
ity, in this work we keep the average pairing interaction 
per site, g = (g A + g B )/2, fixed unless otherwise stated, 
and plot T c vs. the difference S = (g A — g B ). We see 
that the numerical discrepancies between the exact eval- 
uation and the window DOS approximation are in fact 
quite large, due to the presence of the van Hove singu- 
larity at the Fermi level in the simple 2D band, together 
with the exponential dependence of T c on the density of 
states, but that the approximation captures the quali- 
tative dependence of T c on inhomogeneity. It is clear 
that, as in other mean field models of inhomogencous 
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pairing, the inhomogeneity enhances the critical temper- 
ature. We return in Sec. IV to the question of when this 
result breaks down. 

If, instead of holding the average interaction g fixed, 
we fix the attractive interaction on one sublattice g A and 
progressively decrease g B , we find a rapidly decreasing 
T c (see insert to Fig. 3). Depending on the magnitude 
of g B , we either see only an exponential decrease of T c 
or, for very small g B , an instability to the metallic state 
at a critical value determined by the change of sign of 
the effective interaction. This determines the solid line 
in the phase boundary corresponding to the transition to 
the normal state shown in Fig. 1. 

III. QUASIPARTICLE STATES AND 
THERMODYNAMICS 




T/t T/t 




5/t 8/t 



A. Order parameters 

We now consider the order parameters A Q which de- 
velop below T c on the two sublattices. It turns out that 
there are two different cases, which we discuss separately, 
corresponding to coupling constants on the two sublat- 
tices of identical or opposite sign. In general, we can 
write using Eqs. (15), 



A^_ A^_ v 



tanh 



f3E 2 



— tanh 



{PEi 



V 2 



(18) 



We first assume both coupling constants positive, < 
9b < 9 A- If both gaps are positive, then E\,Ei are 
positive. Therefore in limit T — ► the right hand side of 
Eq. (18) will vanish, which gives the solution: 



A B 



9 L 



(19) 



which is exact for the case of both g's attractive. Now 
if we employ the window density of states given by Eq. 
(16), we obtain 



St 



9 A + 9 B sinh I 



a A +g i 



(20) 



For g A = g B , within weak coupling approximation 
(g A ' B «8t), this gives A/fcsTc = e 7 /tt, which is the 
standard BCS result. 

The exact solution for the temperature dependence of the 
gap is shown for the case of weak inhomogeneity in Fig. 
4(a). 

Now we consider the other case when one of the cou- 
pling constants is negative, 3b < < gA- From (15) 
we see that A B is forced to be negative. The variation 
of gap as a function of temperature found numerically 



FIG. 4: (a) Sublattice order parameters A A ' B vs. T/t for 
weak inhomogeneous pairing interaction S = (g A ~g B ) = 0.4t; 

(b) strong inhomogeneity cases, with 8 — 2At (note g B = 
-0.2t). (c) Ratio A Q (T = 0)/T c vs. 5 for a = A, B (solid 
line: full numerical evaluation; dashed-dotted line: wDOS); 
(d) Comparison of full numerical evaluation of A a (T = 0) 
(solid line) with wDOS approximation (dashed-dotted line). 
For all plots, the average coupling is fixed to g — 1. Note in 

(c) ,(d) that the gapless phase discussed in Sections IIIB and 
IIIC correspond, for the parameters given here, to 2 < 5/t. 



and displayed in Fig. 4(b) is somewhat unusual. As 
the temperature increases from zero, the gap increases 
in magnitude and then after attaining a maximum value 
begin to decrease. 

In the case of opposite sign coupling constants, relation 
(19) is no longer valid. In this case one can obtain an ap- 
proximate analytical expression for |A S |/A using the 
wDOS (16), but it is quantitatively inaccurate. A more 
accurate approximation is obtained by replacing the ex- 
act contour of zero quasiparticle energy with a square of 

side 2V2 cos -1 | ^ A "J A - ) . This turns out to be equiv- 



alent to a renormalized wDOS approximation where the 
bare bandwidth 8t is replaced by 47rf. We find 



l A I 




(21) 



In Figs. 4 (c) and (d) we illustrate the utility of 
the wDOS and renormalized wDOS approximations with 
which analytical results were obtained above. While Fig. 
4(d) shows that the inaccuracies arising from the wDOS 
also influence the value of the T = order parameters 
on the two sublattices, the qualitative tendencies with 
S are well reproduced. Furthermore, the ratios A A ' B /T c 
agree quantitatively with the exact numerical results over 
a wide range of inhomogencities (Fig. 4(c)), as the errors 
introduced in the DOS arc largely cancelled by taking 
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FIG. 5: (a) Gaps A A ' B /(4t) vs. T/(4i) for special case 8 = 2t, 
g=t (g B = 0); (b) Condensate amplitude $ s = (c B c B ) for 
same case as (a); (c) Variation of $ s with inhomogeneity 8/t 
for three different temperatures T/t = 0, 0.11, 0.18i. 



the ratio. 

The normal state region shown in the phase diagram 
of Fig. 1 represents the metallic phase where, within the 
current mean field treatment, the only solution found has 
both order parameters A A ' B = 0. 




n/<2 



B. Condensate amplitudes 

There are several measures of the "strength of su- 
perconductivity" in a material which become distinct 
when the material is inhomogeneous. One is the or- 
der parameter, which we have discussed above. A sec- 
ond is the pair wave function, or condensate amplitude 
$ Q = (c^ k „ a c£ a ), such that A a = g a <f> a . To illustrate 
the differences, it is useful to consider first an interesting 
if artificial special case where the coupling is attractive 
on one sublattice and vanishes on another. We there- 
fore take g A = 2t and g B = 0, and plot both the order 
parameters and condensate amplitudes in Fig. 5. 

Using the wDOS, the condensate on either the A or B 
sublattice for the g B = case can be obtained analyti- 
cally at T = as 

, - . , „ , (q = A 7 B, g B = 0) 

(22) 

In other words, although the pairing interaction on the 
B sublattice is zero, the proximity effect is so strong that 
a uniform condensate fraction is induced over the entire 
system. This degeneracy is broken at finite temperatures, 
where the B condensate is weakened relative to that on 
the A sublattice. Of course, the order parameter on the 
B sublattice is zero at all temperatures due to the ab- 
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FIG. 6: The quasiparticle energy (in units of i) in momentum 
space along k x (on x- axis), with k y — for average cou- 
pling g — t and various values of inhomogeneity 8 = g a — g b . 
(a)Spectrum when g A = g B . Dashed curve represents normal 
state g A = g B = for reference, (b) Spectrum for 8 = 0.4i 
(c) Spectrum for 8 = 2t, at transition to gapless state, (d) 
Spectrum for 8 — 2 At, in gapless state. All spectra are shown 
at temperature T = 0.104i. 



sence of coupling on these sites. These differences are 
illustrated in Fig. 5(a) and (b). It is also interesting 
to study the dependence of the subdominant condensate 
on inhomogeneity. The point where one of the coupling 
constants changes sign is clearly a singular point of the 
theory, such that at T = this condensate is immediately 
depleted as soon as the subdominant coupling (B sublat- 
tice) becomes negative. We now discuss the reasons for 
this abrupt change in behavior. 



C. Dispersion and density of states 

In Eqs. (8), we have given the quasiparticle eigenen- 
crgics for this simple model; we now plot them in Fig. 6 
along k x for various representative cases. Of course the 
full quasiparticle spectrum retains the fourfold symmetry 
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FIG. 7: Local density of states (LDOS) at sites A (red) and 
B (blue), (a) LDOSin homogeneous case, 5 = 0; (b) 8 = OAt, 
g = t; (c) <5 = 2t, g = t; (d) 5 = 2.4i, g = t. The dashed line 
represents the normal state 2D tight binding band. 
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FIG. 8: Specific heat C(T)/T of inhomogeneous system vs. 
T/(At). Solid line: BCS homogeneous superconductor with 
2D tight-binding band. Circles: weak homogeneity 8 = 0.8£. 
Triangles: strong homogeneity 8 = 2.5t. Average coupling 
fixed to g — t. 



of the underlying lattice. Panel a) shows the dispersion 
in both the normal and homogeneous BCS case, respec- 
tively. Adding a small inhomogeneity with two attractive 
couplings leads to a splitting of these states, into four 
distinct eigenenergies (panel (b)), thereby reducing the 
overall spectral gap in the system. When one of the cou- 
plings passes through zero, a level crossing occurs at zero 
energy at the zone face (panel (c)). It is this level cross- 
ing which is responsible for the singular behavior in some 
observable quantities as one of the couplings becomes re- 
pulsive. It is accompanied by a filling of the gap by the 
states near the zone face. When the two couplings gA 
and qb have opposite signs, the eigenvalues change sign 
over the Brillouin zone. For example, in panel (d) of Fig. 
6, Ei is always positive, but E\ is positive except in a 
narrow range —yj A A \A B \ < £ < X /A A \A B \. 

Gaplcss behavior in s-wave superconductors with mag- 
netic impurities 14 can be associated with individual low- 
energy quasiparticle resonances, which arise near each 
magnetic impurity and then overlap to form an impurity 
band. The band appears as a residual density of states 
for sufficiently high impurity concentration Similarly, the 
gapless superconductivity arising due to microscopically 
inhomogeneous sign-changing order parameters, which is 
found in the present paper, is reminiscent of respective 
low-energy Andreev states around so-called off-diagonal 
or phase impurities 11,12 . One might imagine an ordi- 
nary s-wave superconductor into which a dilute lattice of 
phase impurities is introduced. As the lattice constant 
of such bound states is decreased, the bound state wave 
functions will interfere and the bound state energies will 
split, leading to a kind of gaplcss superconductivity re- 
lated to that we find here for the dense case. 



To observe this directly and trace the spectral features 
back to their respective sublattices, we construct the local 
density of states on A and B sites, 

P a {u) = £ [|<| 2 <5( W - En) + \v?\ 2 5(to + E n )] . (23) 

Fig. 7 exhibits the LDOS on the two sublattices. In 
panels (a) and (b), when both couplings are positive, we 
clearly see a full spectral gap. But when one of the cou- 
plings tends to zero, as in (c), a sharp peak develops at 
the Fermi level on the associated sublattice. The system 
remains gapless as the coupling constant on this sub- 
lattice is made negative, as in (d). It should be noted, 
however, that the subgap states are not located at zero 
energy, as seen in Fig. 6d, where the dispersion is seen to 
flatten at the zone edge away from zero. This is reflected 
in a narrowly split resonance in the LDOS, as seen in Fig. 
7. 



D. Specific heat 

To exhibit the consequences of the sublattice pairing 
inhomogeneity for thermodynamics, we calculate the spe- 
cific heat from the temperature derivative of the entropy, 



Here S for the free quasiparticles is given by 

S= (/n log (/n) + (1 -/«) log (1 -/«))), (25) 

k,n=l,2 

where /„ = f(E n (k)), and the sum runs only over posi- 
tive eigenvalues. In Figure 8, we plot these results for 
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a few examples. For a superconductor with a small 
amount of inhomogeneity-again keeping the average cou- 
pling fixed-it is clear that the critical temperature in- 
creases relative to the BCS case, the low-T temperature 
dependence remains exponential. On the other hand, 
when one of the coupling constants goes negative, the fig- 
ure shows the existence of the residual density of states 
via a linear term in the specific heat at the lowest tem- 
peratures. 

It is clear that the fractional jump in the specific heat, 
reflecting the condensation energy, decreases with in- 
creasing inhomogeneity, and is approximately constant 
for the gapless regime. Note that this ratio for the ho- 
mogeneous case is not the conventional BCS value of 1.43, 
due to the non-parabolic band used here. 



IV. SUPERFLUID DENSITY 

As Aryanpour et al. have emphasized 6 , models of inho- 
mogeneous pairing can lead to a variety of ground states, 
including insulating ones. To show that the states we 
discuss here arc indeed superconducting, we use the cri- 
teria developed by Zhang et. al. 13 for lattice systems. 
We calculate the supcrfluid density n s , which is the sum 
of the diamagnetic response of the system (kinetic energy 
density) and paramagnetic response (current-current cor- 
relation function) 

n s /m = (~k x ) - A(w = 0, q ^ 0). (26) 

The current-current response is defined as 

A xx (q,iiJ m ) = dTe^ T (jF(q,T)jF(-q,0)) (27) 
Jo 

where the current at i th site is given by 

jf« = it£ ( c l +x ,a^ - 4, CT Ci+*, CT ) • (28) 
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FIG. 9: Superfluid density n a vs. inhomogeneity 8 with fixed 
average coupling g = t at different temperatures T — (solid) 
O.lt (dashed), and 0.18i (dotted line). The range of inhomo- 
geneities 2 < 5/t corresponds to gapless superconductivity. 



The other term we need to evaluate is the expectation 
value of the lattice kinetic energy density operator, which 
in the homogeneous BCS case is directly proportional to 
the superfluid weight n B /m = (—k x ) at T = 0. The 
kinetic energy operator is defined as, 



k x (i) — [ C l+x,cr C h<r + Ci !(T Ci+x,o- 

cr 

In Fourier space, it can be written as: 
k x (i) = -t £ (<t c | e ^^)- 



(30) 



k\ ,/c2 ,cr 



k 2 y \ i£2s. 



(31) 

After performing the canonical transformation and tak- 
ing the expectation values we find 



h.c. 



After evaluating the expectation value using the Bo- 
goliubov operators in (7), we find the following relatively 
simple expression for the analytic continuation of the 
static homogeneous response: 



A ra (w = 0,q^ 0) 



4t sin(— cos(-^=) 
V2 V2 . 



'(29) 



f fjE,) - f(E 2 ) 
V Ek-Ek 



If the Ei(k) do not change sign over the Brillouin zone, 
as is the case for the gapped phase gA , .9 b > 0, it is 
clear that this expression vanishes as T — > as in the 
clean BCS case. This is no longer the case in the gapless 
regime gs < 0, where a finite value of A, corresponding 
to a residual density of quasiparticles, is found at zero 
temperature. 



(Mi)> = -2*£cos(^| 

X E { [< U n 

n 

+ 



B* i— 5= , A* B — i— f= 

e + v £ v"e V2 



f(E n )+ 

(1 -/(£„))} (32) 

which is, remarkably, the same for sites A or B, corre- 
sponding to the uniform proximity-induced condensate 
amplitude found in section IIIB. In terms of the eigen- 
vectors and eigenvalues given by Eqs. (9), (10), (11) and 
(12), we can write 
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FIG. 10: The upper graph shows the superfluid density n s as a 
function of temperature for various inhomogeneity parameters 
8/t= (black), 1.2 (red), 2 (green) and g = t (magenta). The 
lower panel shows the current current correlation A XX (T). 



FIG. 11: Mean field critical temperature T c (dashed-dotted 
line) and phase ordering temperature Tg (dashed line) vs. in- 
homogeneity S for fixed average interaction g = t. The solid 
line is the minimum of the two curves at any S. 



For the case when g A , g B > at T = 0, we can simplify 
this expression as 

(-fc a (i))=2^2 ^ (34) 

S y/ + +4% 

which is manifestly positive; hence, the system is a su- 
perconductor and displays a conventional Meissner effect 
at T = 0. The expression also shows explicitly that the 
superfluid density on each site corresponds to the average 
superconducting order parameter over the system. 

Fig. 9 show how the superfluid weight (n s ) changes as 
the inhomogeneity is increased. At T = 0, n s for this 
model is a constant and equal to the value for the ho- 
mogeneous system, whenever the system is fully gapped, 
since the average gap remains the same. In the same sys- 
tem, recall, the transition temperature increases mono- 
tonically with inhomogeneity. As one increases the inho- 
mogeneity further with fixed average coupling, there is a 
critical value 8 C for which the system enters the gapless 
phase, at which the superfluid density drops abruptly. 
The temperature dependence of such a case is exhibited, 
along with others, in Fig. 10, and we see that this dis- 
continuity corresponds to the creation of a finite residual 
normal fluid of uncondensed quasiparticlcs. Exactly at 8 C 
(= It in Fig. 10), this behavior is marginal, but for any 
8 > 2t in the Figure, the residual normal fluid density 
h-xx{T = 0) is nonzero. 

If g is held fixed and inhomogeneity 8 increased further, 
there is no mean-field transition to the metallic state, as 
indicated e.g. by the dashed-dotted line in Fig. 1 which 
we have investigated in this paper. Instead, T c increases 
monotonically while the superfluid density begins to fall 
at T = when the system enters the gapless phase at 8 C . 
Eventually, the decrease in n s must lead to large phase 



fluctuations and the breakdown of the mean field theory 
described here. To estimate when this occurs, we follow 
Emery and Kivclson 17 and define a characteristic phase 
ordering temperature 

T e = H 2 n s (0)e/4m*, (35) 

which will now decrease in the gapless phase, following 
the superfluid density. For quasi two-dimensional super- 
conductors, the length scale I is a larger of the average 
spacing between layers d or y^Cx, where £j_ is the su- 
perconducting coherence length perpendicular to the lay- 
ers. Although a determination of I is beyond our simple 
model, it is qualitatively irrelevant in the presence of sig- 
nificant abrupt drop of the superfluid density. Taking 
simply t ~ a, which is within a factor of 2-3 of both d 
and v^C-L m the cuprates, we may now obtain a crude 
measure of the optimal inhomogeneity needed to max- 
imize the true ordering temperature in this model. In 
Fig. 11, we now plot both the mean field T c and the pu- 
tative phase ordering temperature as a function of inho- 
mogeneity 8. We see that the two curves cross very close 
to the phase boundary 8 C (at 8 = 2t in the Figure), due to 
the steep drop in n s there. It therefore appears that the 
optimal inhomogeneity within the current scheme occurs 
for a checkerboard-like pairing interaction with attractive 
interaction on one sublattice and zero interaction on the 
other. 



V. CONCLUSIONS 

We have presented and solved a simple mean-field 
model of a superconductor with inhomogencous pairing 
varying on an atomic scale. The model simply takes a dif- 
ferent value for pairing on each of two interpenetrating 
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sublatticcs. The inhomogeneity can be increased while 
keeping the average pairing in the system fixed, thus al- 
lowing the effect of inhomogeneity on various properties 
to be discussed. Many results for this model were pre- 
sented in analytic form. As pointed out by earlier work- 
ers, we find that increasing pairing inhomogeneity en- 
hances T c . We have presented analytical results reflect- 
ing this effect, and calculated the spectrum and thermo- 
dynamic properties as well. The enhancement of T c is 
accompanied by a closing of the spectral gap and eventu- 
ally by the appearance of gapless phase in the quasipar- 
ticle spectrum of the "s-wave" superconductor. While 
these results have been obtained at the mean field level, 
because the length scale of the variation of the pairing 
interaction is so short we do not expect fluctuations to 
change them significantly 3 , at least until the supcrfluid 
density is substantially suppressed. 

The last caveat does not seem to be particularly re- 
strictive, however. The enhanced T c found in the inhomo- 
geneously paired system is accompanied by a reduction 
in condensation energy, but the system maintains a con- 
stant T = supcrfluid density until a critical value of the 



inhomogeneity, when a level crossing of two quasiparticlc 
energies forces the system into the gapless phase. Until 
this point, the condensate at T = is spatially homo- 
geneous for any value of the interaction inhomogeneity, 
even when the interaction on a given sublattice is zero; 
this reflects the strong proximity coupling between the 
sublattices. As inhomogeneity increases further, we have 
argued that the decrease in the supcrfluid density leads 
fairly rapidly to a decrease in the phase ordering temper- 
ature; we find, therefore, that the optimal inhomogeneity 
is reached close to the checkerboard pattern of interac- 
tions, with attractive pairing on one sublattice and zero 
on the other. 

Our results suggest that in superconducting systems 
with short coherence length, a modulated pairing inter- 
action at the atomic scale may provide a route to high 
temperature superconductivity. 
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